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The Influence of Different Turbulence Models on the Diffusion 
Coefficients of Energetic Particles 

M. Hussein^, R. C. Tautz^, and A. Shalchi^ 

Abstract. We explore the influence of turbulence on the transport of energetic parti¬ 
cles by using test-particle simulations. We compute parallel and perpendicular diffusion 
coefficients for two-component turbulence, isotropic turbulence, a model based on Goldreich- 
Sridhar scaling, noisy reduced magnetohydrodynamic turbulence, and a noisy slab model. 

We show that diffusion coefficients have a similar rigidity dependence regardless which 
turbulence model is used and, thus, we conclude that the influence of turbulence on par¬ 
ticle transport is not as strong as originally thought. Only fundamental quantities such 
as particle rigidity and the Kubo number are relevant. In the current paper we also con- 
hrm the unified nonlinear transport theory for noisy slab turbulence. To double-check 
the validity and accuracy of our numerical results, we use a second test-particle code. 

We show that both codes provide very similar results confirming the validity of our con¬ 
clusions. 


1. Introduction 

There are two fundamental problems in plasma physics 
and astrophysics, namely the understanding of magnetic tur¬ 
bulence and the propagation of energetic particles such as 
cosmic rays through a plasma. It is well-known and ac¬ 
cepted that these two problems are linked to each other. In 
analytical treatments of particle transport, for instance, a 
crucial input is the so-called magnetic correlation tensor de¬ 
scribing the turbulence. Energetic particles are assumed to 
move diffusively and, thus, a diffusive transport equation is 
usually used to model their motion (see, e.g., Schlickeiser 
[2002] for a review). Fundamental quantities in that trans¬ 
port equation are the spatial diffusion coefficients, besides 
other transport parameters such as the coefficient of stochas¬ 
tic acceleration. We would like to emphasize that in recent 
years non-diffusive transport has been discussed as well (see, 
e.g., Zimbardo et al. [2006], Pommois et al. [2007], Shalchi 
and Kourakis [2007], Tautz and Shalchi [2010], Zimbardo et 
al. [2012]). 

The turbulent magnetic fields SB, which control the 
transport of energetic particles, are superposed by a mean 
magnetic field Bq. In the general case one has to distin¬ 
guish between diffusion along and across that mean field. In 
the former case the important quantity is the parallel diffu¬ 
sion coefficient ac[| and in the latter case the perpendicular 
diffusion coefficient k±. Together with the so-called drift 
coefficient ka, they form the diffusion tensor. 

The knowledge of diffusion parameters is essential for dif¬ 
ferent applications such as the description of diffusive ac¬ 
celeration at interplanetary shocks (see, e.g., Zank et al. 
[2004], Bosch and Shalchi [2010], Li et al. [2012], Wang 
et al. [2012]). In the context of interstellar shocks, such as 
supernova remnant shocks, usually simple diffusion models 
are employed but it was shown recently that the form of the 
diffusion coefficient can have a strong impact on the cosmic 
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ray spectrum (see Ferrand et al. [2014]). Analytical forms 
of the diffusion coefficients are also required for studies of 
solar modulation and space weather (see, e.g., Alania et al. 
[2013], Engelbrecht and Burger [2013], Manuel et al. [2014], 
Potgieter et al. [2014], Zhao et al. [2014], Engelbrecht and 
Burger [2015]). To understand the propagation of cosmic 
rays in different astrophysical environments, the knowledge 
of transport parameters is also required (see, e.g., Shalchi 
and Biisching [2010], Thornburg and Drury [2014]). 

As mentioned above, a critical quantity entering diffusion 
theories is the magnetic correlation tensor describing the 
turbulence. Therefore, one could conclude that a detailed 
understanding of turbulence is crucial in the theory of en¬ 
ergetic particles. However, there are also some indications 
obtained from analytical theory that transport is universal 
meaning that the details of turbulence are not important 
if diffusion coefficients are calculated. Based on the uni¬ 
fied nonlinear transport (UNLT) theory (see Shalchi [2010]) 
it was shown in Shalchi [2014] that diffusion coefficients of 
magnetic field lines and perpendicular diffusion coefficients 
of energetic particles only depend on fundamental length 
scales of the turbulence and the magnetic field ratio SB/Bq. 
This statement was strengthened by the work of Shalchi 
[2015] where it was shown that the field line diffusion coef¬ 
ficient is only controlled by the Kubo number and the per¬ 
pendicular diffusion coefficient of energetic particles is only 
controlled by the Kubo number and the parallel mean free 
path. The details of the spectrum such as spectral shape or 
anisotropy are not important as long as the aforementioned 
scales are finite (see, e.g., Matthaeus et al. [2007], Shalchi 
[2014], and Shalchi [2015] for details). 

In the current paper we are trying to achieve the follow¬ 
ing: 

1. We perform test-particle simulations for five differ¬ 
ent turbulence models, namely the two-component model, 
isotropic turbulence, Goldreich-Sridhar turbulence, the 
Noisy Reduced MagnetoHydroDynamic (NRMHD) model, 
and a noisy slab model. All models are briefly presented in 
Section 2 of the current paper. We compute numerically the 
parallel and perpendicular diffusion coefficients and compare 
them with each other. We aim to investigate the influence 
of the turbulence onto the different transport parameters. 

2. In Shalchi [2015] a new perpendicular transport regime 
was discovered for the case of short parallel mean free paths 
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and small Kubo numbers. This new regime is very simi¬ 
lar compared to the scaling originally obtained by Rechester 
and Rosenbluth [1978]. In the current paper we simulate for 
the hrst time particle transport in noisy slab turbulence and 
test numerically the aforementioned predictions. 

3. We employ a second test-particle code, namely the so- 
called PADIAN code developed by Tautz [2010] to compute 
the different diffusion coefficients independently. This will 
allow us to check the validity of our simulations and to un¬ 
derstand how accurate and reliable our findings are. 

The remainder of this paper is organized as follows. In 
Section 2 we discuss the different turbulence models em¬ 
ployed in the current paper. The test-particle code which is 
used is discussed in Section 3. In Section 4 we show the dif¬ 
fusion coefficients along and across the mean magnetic field 
for the different turbulence models. In Section 5 we use the 
PADIAN code to check the validity of the results presented 
in Section 4 and we end with a short summary and some 
conclusions in Section 6. 


\i l,m — x,y and Piz = P^rn = Pzz = 0. In this partic¬ 
ular model the magnetic field vector as well as the spatial 
dependence are two-dimensional. Above we have used the 
spectrum of the two-dimensional modes g^^{k±). 

A model which is frequently used to approximate solar 
wind turbulence is the so-called two-component model in 
which we have 

pc^m, ^ psl^a, ^ p2D (g) ^ (4) 

This model is supported by solar wind observations (see, 
e.g., Matthaeus et al. [1990], Osman and Horbury [2009a], 
Osman and Horbury [2009b], Turner et al. [2012]), numer¬ 
ical simulations (see Oughton et al. [1994], Matthaeus et 
al. [1996], Shaikh and Zank [2007]), as well as analytical 
work (see Zank and Matthaeus [1993]). To complement the 
two-component model we have to specify the two spectra 
gsiab^k^^'j g'^^(^k±), respectively. For the former spec¬ 
trum we use the form proposed by Bieber et al. [1994] 


2. Synthetic Models of Turbulence 

In analytical diffusion theories and also in test-particle 
simulations, an important quantity is the so-called magnetic 
correlation tensor defined as 

P,r^{k) ^{5Bl{k)SB*rr,{k)). (1) 

Here we have used the ensemble average operator (...). In 
the general case, the latter tensor can also depend on time 
(see, e.g., Schliekeiser [2002] and Shalchi [2009] for reviews). 
In the current paper we use the magnetostatic approxima¬ 
tion. In the theory of particle transport, the latter approxi¬ 
mation can be justihed by considering only particles moving 
much faster than the Alfven speed. Such particle energies 
and rigidities are considered in the present paper. 

It is still unclear what the exact form of the correlation 
tensor really is in different physical scenarios. One would 
expect that turbulence in the solar system is different com¬ 
pared to turbulence in the interstellar space (see, e.g., Hu- 
nana and Zank [2010] or Shalchi et al. [2010] for a more de¬ 
tailed discussion of this matter). In the following we review 
hve different models which were discussed in the literature 
before. We do not judge which of these models is more re¬ 
alistic nor do we claim that our list of turbulence models is 
complete. 

2.1. The two-component model 

In very early treatments of particle transport, simple 
models for the turbulence have been employed. Jokipii 
[1966], for instance, used a so-called slab model in which the 
magnetic correlation tensor has by dehnition the following 
form 


E*l“^(fc)=/'“^(fc||)^5*^, (2) 

with l,m = X, y. Here we have used the Kronecker delta 
Sim and the Dirae delta 5{k±), respectively. The other ten¬ 
sor components are zero due to the solenoidal constraint. 
Furthermore, we have used the spectrum of the slab modes 
^sia6(^ll) have used the wave vector components 

along and across the mean magnetic field fc|| and k^, re¬ 
spectively. The slab model is basically a one-dimensional 
model in which the turbulent magnetic field depends only 
on the coordinate along the mean field. 

Another model with reduced dimensionality is the so- 
called two-dimensional model where we have by definition 




^SBl,l,lslab 


[l -|- (fc||Zs 


s/2 ’ 


( 5 ) 


where we have used the slab bendover scale I slab, the mag¬ 
netic held strength of the slab modes SBsiab, and the inertial 
range spectral index s. The normalization function C{s) 
is given below. For the spectrum of the two-dimensional 
modes we employ the one proposed by Shalchi and Wein- 
horst [2009] 


2D N ‘2D{s, q) r.j^2 , 

g {k±) = - SB 2 DI 2 D 


{kj2 


[l + (fc±Z2o)2](^+«)/"’ 


( 6 ) 


where we have used the two-dimensional bendover scale hn, 
the magnetic held strength of the two-dimensional modes 
SB 2 D, and the energy range spectral index q. The normal¬ 
ization function is given by 


D{s,q) 






( 7 ) 


where we have used the Gamma function r( 2 ). The spec¬ 
trum (6) is correctly normalized if s > 1 and q > —1 are 
satished. The normalization function of the slab spectrum 
is a special case of the function D{s,q). The two functions 
are related to each other via C{s) = D{s,q = 0). To en¬ 
sure that the ultra-scale and other fundamental turbulence 
scales are hnite, we only consider cases with q > 1 (see, e.g., 
Matthaeus et al. [2007] and Shalchi [2014]). 


2.2. Isotropic turbulence 

Before a more detailed understanding of turbulence be¬ 
came available, scientists used two simple models, namely 
the slab model discussed above and the isotropic model (see, 
e.g., Fisk et al. [1974] and Bieber et al. [1988]). In the latter 
case the turbulent magnetic held itself depends on all three 
coordinates of space but there is no preferred direction. In 
this case the correlation tensor has the form 


P ,iso 
Im 


g"‘‘°{k) (^ kikm\ 

87 rfc 2 V j ’ 


( 8 ) 


where we have used the isotropic spectrum g’'‘‘°{k) which is 
defined so that 


d 2D/J*\ 2D n 


^Im 


kl krri 


( 3 ) 


dk g'-‘°{k) = SB'^ 


k_ 


(9) 



HUSSEIN ET AL.: DIFEUSION COEFEICIENTS OF ENERGETIC PARTICLES 


X - 3 


For the latter spectrum we employ the form 

( 10 , 

corresponding to the spectrum used above for two- 
dimensional turbulence. In the current paper we set q ~ 3 
if the isotropic model is considered to ensure again hnite 
turbulence scales. The parameter Iq denotes the bendover 
scale. 

2.3. Goldreich-Sridhar turbulence 


the model used in Ruffolo and Matthaeus [2013]. The func¬ 
tion g^^{k±) is the spectrum of the two-dimensioual modes 
as used above. It is obvious from the definition (13) that 
the NRMHD model can be understood as a broadened two- 
dimensional model. We can easily recover the pure two- 
dimensional model discussed above by consideriug the limit 
K^Q. 

In the following we use the same spectrum which was used 
before by Ruffolo and Matthaeus [2013] as well as by Shalchi 
and Hussein [2014] in the context of NRMHD turbulence, 
namely a spectrum with q = 3. In this particular case Eq. 
(6) becomes 


The Goldreich-Sridhar model (see Goldreich and Sridhar 
[1995]) does uot predict the exact form of the tensor of 
Alfvenic perturbations, it only states that most energy will 
be in the state of critical balance. Combining the critical 
balance condition with the Kolmogorov cascading, one gets 
the relation between fey and k±, namely, |fc||| ~ k^^. 

The caveat here is that the wavenumbers are defined in 
the global magnetic field frame of reference, this may be mis¬ 
leading. The critical balance is the condition satisfied only 
in the frame related to the local magnetic field. In what 
follows we use the parallel and perpendicular wavenumbers 
but understand them in terms of wavelets, which allow the 
choice of the local magnetic field. More details concerning 
this matter can be found in Lazarian and Vishniac [1999] 
as well as in Kowal and Lazarian [2010]. 

Cho et al. [2002] have proposed a specific form of the 
magnetic correlation tensor which is in agreement with the 
relations discussed above. Shalchi [2013] has suggested to 
use the following generalization of that form 

( 11 ) 


Z'' (fcx) = 


(fexC 


[l -|- {kMi 


n217/3- 


(14) 


The model described here was already used in transport the¬ 
ory to compute field line diffusion coefficients (see Ruffolo 
and Matthaeus [2013] and Snodin et al. [2013]) and perpen¬ 
dicular diffusion coefficients of energetic particle (see Shalchi 
and Hussein [2014]). An interesting aspect of the NRMHD 
model is the fact that it contains two characteristic length 
scales, namely the perpendicular scale l± and the parallel 
scale Z|| = K~^. It was shown in Shalchi and Hussein [2014] 
and Shalchi [2015] that the scale ratio l\\/l± has a strong 
influence on the perpendicular diffusion coefficient. 


2.5. A noisy slab turbulence 

Above we have used the NRMHD model which can be 
understood as a broadened two-dimensional model. We can 
combine the same idea with the slab model which is done in 
the current paragraph. We define the noisy slab model via 
(see Shalchi [2015]) 


with 




D{s,q)_iS^^2 


2tt 


{kjf 


\l + {k±l) 


21 (s+'j)/2 


( 12 ) 


which is only valid for s = 5/3. Compared to Gho et al. 
[2002], the latter spectrum takes into account a more general 
Isehavior at large scales corresponding to the energy range. 
This spectrum and similar analytical forms were used before 
in the theory of energetic particle transport (see, e.g., Chan- 
dran [2000], Cho et al. [2002], Shalchi et al. [2010], Sun and 
Jokipii [2011], and Shalchi [2013]). In Shalchi [2014] funda¬ 
mental length scales of turbulence were computed for this 
type of spectrum. There it was shown that the ultra-scale is 
only finite if g > 1. Therefore, we set g = 2 if the Goldreich- 
Sridhar model is considered. 


2.4. The NRMHD model 

Above we have discussed the two-dimensional turbulence 
model. One could argue, that this model is a singular 
model because there is no variation of the turbulent field 
along the z-axis. Ruffolo and Matthaeus [2013] proposed a 
so-called Noisy Reduced Magneto Hydro Dynamic (NRMHD) 
turbulence model in which the magnetic correlation tensor 
(1) has the following form 

pr^^(fc) = Q%^e(i^-|fcii|)(5*^-^) (13) 

where we have employed the Heaviside step function 0{x) 
which is defined so that 0{x > 0) = 1 and 0{x < 0) = 0. 
Therefore, Q{K — |fc|||) = 0 for |fc||| > A in agreement with 


9/ 

P,//(fc) = |ig*'“^(fc|l)e(l-fcxZx) 





(15) 


where l± is a characteristic length scale for the decorrela¬ 
tion across the mean magnetic field. In the case that one 
or two indexes are equal to z, the element of the correlation 
tensor is assumed to be zero as in the pure slab model de¬ 
fined above. If there is no broadening, corresponding to the 
case Zx — >■ CX3, the noisy slab model corresponds to the usual 
slab model. In Eq. (15) we have used again the Heaviside 
step function and p^*“*’(fc||) is the usual spectrum of the slab 
modes as it was used above. To study particle diffusion in 
noisy slab turbulence is interesting because Shalchi [2015] 
predicted that the ratio Ax/Ay is much smaller in this case 
compared to other turbulence models. 

3. Test-Particle Simulations 

In the current article, we perform simulations to obtain 
parallel and perpendicular diffusion coefficients numerically. 
In computer simulations three steps have to be performed 
in order to obtain diffusion parameters, which are 

1. A specific turbulence model has to be simulated by em¬ 
ploying the approach described below. Here, we consider the 
slab/2D composite model, isotropic turbulence, Goldreich- 
Sridhar turbulence, the NRMHD model, and a noisy slab 
model. 

2. The Newton-Lorentz equation has to be solved numer¬ 
ically for an ensemble of particles to obtain their orbits. 

3. From these test-particle trajectories, one can obtain 
the diffusion coefficients in the different directions of space. 

This method of generating the turbulence and simulating 
the motion of test-particles was used before (see, e.g., Gi- 
acalone and Jokipii [1994], Michalek and Ostrowski [1996], 
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Reville et al. [2008], Tautz [2010]) and is different compared 
to the grid method used by other authors (see, e.g, Qin et 
al. [2002a] and Qin et al. [2002b]). 

3.1. General remarks 

In order to calculate the turbulent magnetic field at the 
position of the charged particle x, one can use the Fourier 
representation 


SB{x)^ J d^k5B{k)e*^ ". (16) 

In order to benefit from symmetry, it is preferred to evalu¬ 
ate this integral either in spherical or cylindrical coordinates. 
Since we are dealing with a numerical treatment, integrals 
are replaced by sums. The basic idea is to generate ran¬ 
dom magnetic fluctuations by superposing a large number 
of plane waves with different and random polarizations and 
phases. The dimensionality of the used turbulence model 
matters since it determines to how many sums the integral 
will break up. For example, in turbulence models with re¬ 
duced dimensionality, such as slab or two-dimensional mod¬ 
els, and for isotropic turbulence, the integral can be replaced 
by a single sum because only one independent wave vec¬ 
tor component controls the turbulent magnetic field. On 
the other hand, Goldrich-Sridhar, NRMHD, and noisy slab 
models are more complicated. This is due to the fact that 
two wave vector components are relevant, namely fc|| and 
k±. Therefore, an extra sum is required making the simu¬ 
lations more time consuming. In the following we elaborate 
on the technical details for the different turbulence models. 

3.2. Slab, two-dimensional, and isotropic turbulence 

In order to simulate turbulence models with reduced di¬ 
mensionality, we use the method described in Hussein and 
Shalchi [2014a] and Hussein and Shalchi [2014b]. More de¬ 
tails about the numerical approach used in such simulations 
can be found in Tautz and Dosch [2013]. In numerical treat¬ 
ments of the transport, it is convenient to use dimensionless 
quantities instead of the physical quantities used above. If 
the slab model is considered, for instance, physical quanti¬ 
ties are the parallel component of the wave vector fc|| and the 
parallel particle position a. In the test-particle code those 
quantities are replaced by fc|| —>• k\\lsiab and 2 zlhiab, re¬ 
spectively. More details concerning other quantities such as 
the magnetic rigidity are discussed below. 

In numerical work the turbulent magnetic field vector is 
calculated via 


N 

SB (x) = V2SB ^ A{kn)in cos [k„ ■ x + l3„] . (17) 

n=l 

The parameter N corresponds to the number of simulated 
wave modes and the quantity A{kn) denotes the amplitude 
function (see below). In Eq. (17) we have also used the 
wave vector kn = knkn with the random wave unit vector 

( \/ l - rjl cos ((>„ \ 

\/l-rin sin (j)„ j ■ (18) 


Furthermore, we have used the random phase fdn and the 
polarization vector 


In 


— sin 4>n cos an + rjn cos (j>n sin a 
cos 4>n cos an + Tjn Siu <()n Siu 


■A 


- rjf, sm an 


(19) 


with rjn = cos^n. The angles 9n, (j>n, and an can have a 
specific value or they are random angles depending on the 
simulated turbulence model (see Table 1 of the current paper 
for the used values). 

The amplitude function A{kn) used above depends on the 
spectrum G{kn) via 


A^kn) = G{kn)Ak 



G(fc^)Afc^ 


( 20 ) 


For the spectrum we use a form corresponding to the ana¬ 
lytical models described above, namely 

The parameters q and s are energy and inertial range spec¬ 
tral indexes as described above. For the slab modes we use 
q — 0, for two-dimensional modes q = 2, and for isotropic 
turbulence q = 3. A Kolmogorov [1941] spectrum with 
s = 5/3 is used in all cases. To simulate slab turbulence 
we set the polar angle = 0 corresponding to = 1. For 
two-dimensional turbulence, we set r/n ~ 0 and = 0. For 
isotropic turbulence all angles are randomly generated. The 
used values and the meaning of in the different models 
are summarized in Table 1. The composite model is cre¬ 
ated upon superposing slab and two-dimensional modes via 
Eq. (4). In the current paper we set SB^i^f^/Bg = 0.2 and 
SB 20 /B 0 = 0.8 as suggested by Bieber et al. [1996]. 

Eq. (17) is used to compute the magnetic field at the 
considered position. An alternative approach to compute 
turbulent magnetic fields is based on a Fast Fourier Trans¬ 
form to replace the summation (see, e.g., Decker and Vlahos 
[1986], Decker [1993]). 

3.3. Anisotropic three-dimensional models 

In the current paper we also simulate anisotropic three- 
dimensional models such as the Goldreich-Sridhar model, 
NRMHD turbulence, and the noisy slab model. In such 
cases we replace the single sum in Eq. (17) by a double 
sum. Now the turbulent magnetic field vector is given by 


SB{x) = 


X 


M N 

V2SB ^ ^ A{kn, km)in 

m=l n=l 

cos {kn COS (fnX -f fc„ siu (fnV + kmZ Pn,m) ( 22 ) 


corresponding to cylindrical coordinates where kn represents 
the perpendicular wave number k± and km represents the 
parallel wave number fc||. In Eq. (22) we have used again 
the polarization vector fn given by Eq. (19) but set Tjn ~ 0 
therein. The parameters M and N denote the number of 
wave modes in parallel and perpendicular directions, respec¬ 
tively. In the case of NRMHD and noisy slab turbulence, we 
set On = 0 to ensure that SBz = 0. For Goldrich-Sridhar 
turbulence, however, an is a random angle. 

The function A{kn,km) used in Eq. (22) represents the 
wave amplitude and parameter /3„,m denotes the random 
phase as before. For the amplitude function A{kn,km) we 
employ 


A^{kn, km) 


G{kn ) kn Akm Akn 

E^=lEtlGikn)knAk,Ak, 


(23) 


and the turbulence spectrum G{kn) is in the case of the 
NRMHD model given by 


G{kn) = 


kl 


-1 


(l + fc2)(-+^)/2- 


(24) 
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In the current paper we set q — 3 as originally used in Rujfolo 
and Matthaeus [2013]. Furthermore, we cut off the spectrum 
in the parallel direction by using a maximum wave number 
corresponding to the parameter K used in Eq. (13). 

For the turbulence model based on Goldrich-Sridhar scal¬ 
ing we employ the spectrum 


G{km km) — 


kt 


(1 + fe 2 )(-+ 0/2 


-\km\K 


(25) 


with s = 5/3 and q = 2 as discussed above. 

For the noisy slab model, we use the same spectrum used 
for slab turbulence, but we cut off the spectrum in the per¬ 
pendicular direction by using a maximum wave number cor¬ 
responding to the parameter l\\/l±- 


3.4. Further parameters and accuracy issues 

For all anisotropic three-dimensional models, Akm and 
Akn are the spacings between wave numbers. In our simu¬ 
lations we use a logarithmic spacing in km and so that 


Afc„ 

kn 


exp 


ln(fc 

n,max jkn^min ) 

N- 1 


(26) 


and the same for km- It is important that parallel wave 
numbers are distributed fine enough so that the so-called 
resonance condition is satisfied. The resonance condition 
occurs in quasilinear treatments of the transport and states 
that parallel scattering occurs only if jaR^k^^ — 1. In the lat¬ 
ter condition we have used the unperturbed Larmour radius 
Rl at /i = 0 and the pitch-angle cosine /i. In the simulations 
we have to ensure that a large amount of wave numbers are 
close to the corresponding fc||. 

The size of the box is restricted by the so-called scaling 
condition that ensures that no particles travel beyond the 
maximum size of the system, Lmax = This is ensured 

via the relation MmaxkminRh < 1, which corresponds to 
vtmax < Lmax- lu parallel and perpendicular directions we 
used kmin = 10~®, leading to a relatively large box to ensure 
that finite box size effects do not occur. For the maximal 
wavenumber we used kmax = 10® in both directions. 

The simulations contain further parameters controlling 
the accuracy. We have to specify the number of wave modes 
in the parallel direction N and perpendicular direction M, 
respectively. For the simulations performed for slab, two- 
dimensional, and isotropic turbulence, there is only one 
maximum wave number. In this case we have used N = 512 
which is enough to satisfy the aforementioned conditions. 
For the anisotropic three-dimensional models, the double 
sum makes it more challenging concerning computational 
time. For the NRMHD model, we have used N = 256 and 
M = 32, respectively. We have performed test runs with M 
up to 128 and no significant differences were noticed. For 
the simulations performed for Goldrich-Sridhar and nosiy 
slab turbulence we have used N = 32 and M = 256. 

The procedure described so far can be used to compute 
the turbulent magnetic field at the position of the charged 
particle. To obtain the trajectories of the energetic and 
electrically charged particles, we solve the Newton-Lorentz 
equation numerically for 1000 particles. The latter equation 
has the form 


^=^FxB[f(t)] (27) 

dt c 

where x{t) is the position of the particle in Cartesian coor¬ 
dinates. In Eq. (27) we replace the turbulent magnetic field 
by using either Eq. (17) or (22). It is worth noting that we 


are using the relativistic version of the Lorentz force and so 
the momentum is given by p = 'ymv. The simulations are 
done using a Monte Carlo code where each particle is given a 
random initial position, pitch angle cosine /i, and turbulence 
angles. After injection, particles are traced for a sufficiently 
long time, around tens of thousand of gyro-periods for the 
particle to overcome the ballistic regime and move diffu¬ 
sively. For our numerical integrator we use Runge-Kutta of 
4th order which keeps truncation error relatively small and 
under control. 

We have traced 1000 particles for a maximum running 
time of Tmax = Lltmax = u X 10"^ with integer n £ 1, 2,..., 10 
depending on the particle rigidity. When the particle is less 
energetic, i.e. have lower rigidity, it takes more time to move 
diffusively especially in the perpendicular direction. After 
that we have averaged over all the 1000 realizations we have 
been using. To estimate the error calculated for the diffu¬ 
sion coefficients, we follow the method explained in Tautz 
[2010]. This method takes into account the averaging pro¬ 
cedure over turbulence manifestation and over all particles. 
Within our calculations the error was relatively small given 
the fairly large number of turbulence modes integrated over 
and the number of realizations used. This allowed for dif¬ 
fusion to be fairly consistent and stable so deviations from 
the mean diffusion coefficient where hardly noticed. For 
that reason we don’t include errors when plotting out final 
graphs. 

4. Diffusion Coefficients obtained from 
Test-Particle Simnlations 

In the current section we perform the test-particle simu¬ 
lations described in Section 3. In order to distinguish these 
results from the simulations presented in Section 5, we re¬ 
fer to the code used here as the simulations by Hussein & 
Shalchi. In the following we compute the parallel mean free 
path A|| which is related to the parallel spatial diffusion co¬ 
efficient fv|| via A|| = 3 k[|/u, as well as the perpendicular 
mean free path A^ and the ratio A^/Ay. In Tables 2-6 we 
summarize our numerical results. For all simulations we 
set 5B = Bo and s = 5/3. The values for the energy range 
spectral index q are shown in the corresponding table. Some 
models contain two scales ly and l±. The used values for the 
ratio l±/li\ are also listed in the corresponding table. 

The main aim of the current section is to explore the 
rigidity dependence of the parallel and perpendicular mean 
free path, respectively. However, to ensure that we indeed 
obtain diffusive transport, we also show the (running) diffu¬ 
sion coefficients versus time. In Fig. 1 we have shown the 
parallel diffusion coefficient versus the dimensionless time fit 
for R = 1. Clearly the diffusion parameters become constant 
in time after the well-known initial ballistic regime. There¬ 
fore, we conclude that transport is indeed diffusive for the 
considered cases. Fig. 2 shows the perpendicular diffusion 
coefficient versus time. Again we find diffusive transport 
after the ballistic regime. 

In Fig. 3 we show the parallel mean free paths for the dif¬ 
ferent turbulence models except the noisy slab model which 
is discussed separately (see below). All quantities are nor¬ 
malized with respect to L which stands for the corresponding 
turbulence scale (e.g., L = lo for isotropic turbulence). In 
all cases the parallel diffusion coefficient increases with in¬ 
creasing rigidity. The rigidity dependence is approximately 
the same in each case and only the absolute values are 
slightly different. The only exception is the parallel mean 
free path obtained for the NRMHD model which is clearly 
larger at high rigidities. We would like to point out that 
in this particular turbulence model, there is no turbulence 
for certain parallel wave numbers. Therefore, there is often 
no gyro-resonant interaction between the energetic particles 
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and magnetic fields. This could explain the differences be¬ 
tween the parallel diffusion coefficients shown in Fig. 3. The 
analytical investigation of parallel diffusion in NRMHD tur¬ 
bulence will be subject of future work. 

The perpendicular mean free paths for the different tur¬ 
bulence models are visualized in Fig. 4. Again we obtain 
similar results for all considered models. In all cases the 
perpendicular mean free path increases with rigidity if the 
latter parameter is small. For high rigidities the perpen¬ 
dicular mean free paths become rigidity independent as al¬ 
ready predicted analytically in Shalchi [2014] and Shalchi 
[2015]. In this case the perpendicular diffusion coefficient 
approaches asymptotically the so-called Field Line Random 
Walk (FLRW) limit. We like to emphasize that 5Bz ^ 0 for 
isotropic and Goldreich-Sridhar turbulence whereas = 0 
for the other models. This could cause a difference in the 
perpendicular diffusion coefficient. 

In Fig. 5 we show the ratio A±/Ay for the different tur¬ 
bulence models. For low rigidities the ratio Ax/Ay is ap¬ 
proximately constant as expected (see again Shalchi ]2014] 
and Shalchi [2015]). Rigidity dependence and magnitude of 
the perpendicular mean free path depend only weakly on 
the chosen turbulence model. An exception is the NRMHD 
model where the parallel diffusion coefficient is different 
compared to the other models. A possible explanation for 
that difference is provided above. 

A direct relation between parallel and perpendicular 
mean free paths is predicted by the UNLT theory of Shalchi 
]2010] where the rigidity does not explicitly enter the corre¬ 
sponding integral equation. Therefore, we show the perpen¬ 
dicular mean free path versus the parallel mean free path 
for the different turbulence models in Fig. 6. As predicted 
by the aforementioned theory, we find Ax ~ Ay for short 
parallel mean free path and Ax = const for long parallel 
mean free path. 

So far we did not discuss the noisy slab model. The reason 
is that this model provides very different results if it comes 
to the perpendicular diffusion coefficient. This difference 
was already predicted analytically by the UNLT theory (see 
Shalchi ]2015]). The latter theory states that in the limit 
A||/f|| —>■ 0 and for small Kubo numbers, the ratio of the two 
mean free paths is given by 


^11 


2 ^11 SB^ 


(28) 


for the noisy slab model. In the current paper we set 
l\\/l± = 0.5, SB^IBq = 1, = 1, and s = 5/3 leading 

to C{s = 5/3) ~ 0.12. In this case Eq. (28) provides 
Ax/A|| « 0.9 X 10“^. For the limit Ay/Zy —>■ oo and small 
Kubo numbers we expect to find the quasilinear scaling (see 
again Shalchi [2015]) and the perpendicular mean free path 
is given by 


/^X on , , 2 ou 

17 = bF 


(29) 


for noisy slab turbulence. For the parameter values used in 
the current paper (see above), this becomes Ax/^y ~ 0.57. 
In Figs. 7 and 8 we show these two analytical limits together 
with the simulations. We can see that the analytical results 
are perfectly in agreement with the simulations. According 
to Table 6 we find that the parallel mean free path is sim¬ 
ilar compared to the other turbulence models. The ratio 
Ax/Ay, however, is much smaller for noisy slab turbulence 
as predicted and explained in Shalchi [2015]. 


5. Testing our Results by using the 
RADIAN Code 

To check the validity and accuracy of the results presented 
in the previous section, we perform test-particle simulations 
also by using a different code. In Tautz ]2010] the so-called 


PADIAN code was developed. This code is similar compared 
to the code described above. A major difference between the 
two codes is that for three-dimensional turbulence the PA¬ 
DIAN code still uses a single sum whereas a double sum is 
used in the Hussein & Shalchi code. A single sum can be 
realized by first determining the two-dimensional surface el¬ 
ement of the unit sphere. This is obtained from the distances 
between the randomly chosen angles in Eq. (19). By sort¬ 
ing the combinations into a matrix with increasing values 
in the rows and columns, the distances between neighboring 
numbers—and thus the two-dimensional volume elements— 
can be calculated. While this approach certainly does not 
’’fill” the surface with values, it is nevertheless a viable al¬ 
ternative to the double sum introduced in Eq. (23). 

Again we compute the different diffusion coefficients and 
compare the numerical findings with the results obtained 
above. In Tables 2-5 we summarize the PADIAN results for 
the individual turbulence models. In Fig. 9 we compare the 
parallel and perpendicular mean free paths obtained from 
the two different simulations for two component turbulence. 
Obviously we find an almost perfect agreement confirming 
the validity of our numerical work. In Fig. 10 the same 
comparison is shown but for isotropic turbulence. Again 
the agreement obtained by using the two different numeri¬ 
cal tools is almost perfect. For anisotropic three-dimensional 
models such as Goldreich-Sridhar turbulence or NRMHD 
turbulence, it is more difficult to perform test-particle simu¬ 
lations and, in this case, there is a major technical difference 
between the two codes. As described above, the Hussein & 
Shalchi code uses a double sum whereas the PADIAN code is 
still using a single sum as in the case of reduced dimension¬ 
ality. In Fig. 11 we compare our results for the Goldreich- 
Sridhar model and Fig. 12 for the NRMHD model. We can 
see that the results are very close although not as close as 
before. Still the agreement allows us to conclude that our 
numerical findings are accurate. Furthermore, it seems to be 
possible to perform the simulations by using a single instead 
of a double sum even if the turbulence is three dimensional. 


6. Summary and Conclusion 

We studied the transport of energetic particles interact¬ 
ing with magnetic turbulence. In order to compute diffusion 
parameters describing the transport, one has to specify the 
properties of the magnetic correlation tensor. In the current 
paper we calculated spatial diffusion coefficients for different 
turbulence models, namely 

• The slab/2D composite model, 

• Isotropic turbulence, 

• A model based on Goldreich and Sridhar [1995] scaling, 

• The NRMHD model of Ruffolo and Matthaeus [2013], 

• A noisly slab model. 

For all those models we computed the parallel mean free 
path Ay, the perpendicular mean free path Ax, and the ratio 
Ax/A||. 

We have shown that for all considered turbulence models, 
the diffusion coefficients have a similar rigidity dependence 
and only the absolute values of the diffusion coefficients are 
different. This conclusion is in agreement with the analytical 
findings obtained in Shalchi [2014] and Shalchi ]2015] based 
on the UNLT theory. The latter theory predicts that the 
perpendicular diffusion coefficient is directly proportional 
to the parallel diffusion coefficient for small rigidities and 
becomes rigidity independent for higher rigidities. This is 
exactly what we found numerically in the current paper (see 
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Figs. 4-8). According to our present work, the influence of 
detailed turbulence properties on the rigidity dependence is 
minor. As already proposed in Shalchi [2014] and Shalchi 
[2015], only fundamental properties of turbulence such as 
the length scales and magnetic fields control the diffusion 
coefficients. As argued in Shalchi [2015], the important 
quantity controlling the transport is the Kubo number. If 
the latter number is extreme, the diffusion parameters are 
extreme as well. This is in particular the case for models 
with reduced dimensionality such as the slab or the two- 
dimensional model. Of course, the statements made in the 
current paper do not apply for such turbulence models. 

In the current paper we have only explored particle trans¬ 
port for certain parameter regimes (e.g., we have assumed 
that SB — Bo). Analytical theories (see again Shalchi 
[2015]) predict that perpendicular diffusion does only de¬ 
pends on the parallel diffusion coefficient and the Kubo num¬ 
ber. The latter number is directly proportional to the ratio 
SB/Bo. Therefore, we expect that changing the magnetic 
field ratio will have a strong influence on the magnitudes 
of the two diffusion parameters. However, qualitatively the 
diffusion parameters should be similar for weak turbulence 
(e.g., for SB/Bo = 0.1 or SB/Bo = 0.3). The latter state¬ 
ment is supported by the numerical work presented in Hus¬ 
sein and Shalchi [2014a] where the influence of the magnetic 
field ratio on the transport had been explored. 

To explore the influence of turbulence on the diffusion 
coefficients of energetic particles was also subject of previ¬ 
ous work (see, e.g., Giacalone and Jokipii [1999], Sun and 
Jokipii [2011]). The numerical results presented in this pre¬ 
vious work is similar compared to our findings. However, 
we have added more turbulence models (e.g., the NRMHD 
model) and we have used a more appropriate spectrum at 
large scales to ensure that all fundamental scales of turbu¬ 
lence (e.g., the ultra-scale) are finite. 

For the first time we have explored numerically transport 
of particles in noisy slab turbulence. For the particular tur¬ 
bulence model we found a very small ratio Ax/Ay in perfect 
agreement with the results obtained in Shalchi [2015]. Once 
more the UNLT theory is confirmed by numerical work. 

To double-check our findings and to estimate how accu¬ 
rate our numerical results are, we have employed a second 
test-particle code, namely the so-called RADIAN code orig¬ 
inally developed in Tautz [2010]. We have shown that this 
second code provides similar results. Only small variations 
can be found if the results provided by the two codes are 
compared with each other. Therefore, we conclude that our 
numerical findings are correct and accurate enough to draw 
conclusions concerning the influence of turbulence on the 
transport of energetic particles. 
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Table 1. The values used in the simulations for slab, two- 
dimensional, isotropic, NRMHD, Goldreich-Sridhar, and noisy 
slab turbulence. 


Turbulence model 

Vn 

Ctrl 

4>„ 

Wave numbers 

Energy range spectral index 

Slab 

1 

0 

Random 

kn — Islab^W 

0 

Two-dimensional 

0 

0 

Random 

kn = l2Dk± 

2 

Isotropic 

Random 

Random 

Random 

kn — lok 

3 

NRMHD 

0 

0 

Random 

kn — l_\_kj^,km ~ ^ll^ll 

3 

Goldreich-Sridhar 

0 

Random 

Random 

kn — ^kj^,kni ~ ^^11 

2 

Noisy slab model 

0 

0 

Random 

kn ~ ^J.k_\_,kn^ ~ 
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Table 2. Test particle simulations for slab/2D turbulence. 
For the energy range spectral index of the two-dimensional 
modes we used q = 2. The two scales are assumed to be equal 
ho = Islab and we set = 0-2 and = 0.8 

as suggested by Bieber et al. (1996). Listed are the results 
obtained by using the Hussein Shalchi code (HS) and the 
PADIAN code (P). 


Rl /^slab 

0.01 

0.05 

0.1 

0.5 

1.0 

5.0 

10.0 

/Islab 

/I slab 

0.90 

1.88 

2.7 

8.5 

15.0 

87.0 

255 

0.048 

0.08 

0.11 

0.26 

0.35 

0.43 

0.45 

A_|_//sZab 

0.932 

- 

2.581 

- 

15.0 

- 

266 

0.049 

- 

0.123 

- 

0.384 

- 

0.456 


Table 3. Test particle simulations for isotropic turbulence. 
For the energy range spectral index we used g = 3 as ex¬ 
plained in the text. Listed are the results obtained by using 
the Hussein & Shalchi code (HS) and the PADIAN code (P). 


Rl/ lo 

0.01 

0.05 

0.1 

0.5 

1.0 

5.0 

10.0 


0.51 

1.05 

1.35 

3.55 

7.1 

82.0 

330 

xr/io 

0.018 

0.04 

0.06 

0.16 

0.235 

0.305 

0.32 

xf/I^ 

0.52 

- 

1.23 

- 

6.58 

- 

325 

\l/lo 

0.018 

- 

0.059 

- 

0.257 

- 

0.319 


Table 4. Test particle simulations for Goldreich-Sridhar tur¬ 
bulence. For the energy range spectral index we used q = 2. 
Listed are the results obtained by using the Hussein Sz Shalchi 
code (HS) and the PADIAN code (P). 


Rl/1 

0.01 

0.05 

0.1 

0.5 

1.0 

5.0 

10.0 

Xf^/l 

0.62 

1.01 

1.34 

3.3 

9.25 

98.0 

410.0 

0.0215 

0.052 

0.0805 

0.255 

0.365 

0.40 

0.45 

xl/l 

0.800 

- 

1.478 

- 

7.186 

- 

386 

0.0251 

- 

0.0715 

- 

0.2711 

- 

0.348 


Table 5. Test particle simulations for NRMHD turbulence. 
For the energy range spectral index we used q = 3 and we 
assumed that /y = l±. Listed are the results obtained by 
using the Hussein Sz Shalchi code (HS) and the PADIAN code 

(p). 


Rl/1± 

0.01 

0.05 

0.1 

1.0 

5.0 

10 

A||//x 

1.05 

2.1 

3.2 

31 

700 

1875 

X±/li_ 

0.04 

0.063 

0.088 

0.195 

0.25 

0.25 

A||//x 

1.279 

- 

5.002 

56.47 

- 

2843 

AxAx 

0.0257 

- 

0.0948 

0.4255 

- 

0.4302 


Table 6. The Hussein Sz Shalchi simulations for noisy slab 
turbulence. For the scale ratio we have assumed = 0.5 

and the used spectrum corresponds to the one used for the 
standard slab model. 


Rl/1\\ 

0.001 

0.01 

0.1 

0.316 

1.0 

3.16 

10.0 

31.6 

100 

A||/^|| 

0.3 

0.66 

1.35 

2.15 

4.15 

12.4 

77 

730 

7700 

X±/l\\ 

0.002 

0.005 

0.0085 

0.012 

0.032 

0.181 

0.40 

0.51 

0.59 
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Figure 1. The running parallel mean free paths from 
the Hussein & Shalchi simulations. Shown are the 
results obtained for slab/2D turbulence (dotted line), 
isotropic turbulence (solid line), Goldreich-Sridhar tur¬ 
bulence (dashed line), and NRMHD turbulence (dash- 
dotted line). 



Figure 2. The running perpendicular mean free paths 
from the Hussein & Shalchi simulations. Shown are the 
results obtained for slab/2D turbulence (dotted line), 
isotropic turbulence (solid line), Goldreich-Sridhar tur¬ 
bulence (dashed line), and NRMHD turbulence (dash- 
dotted line). 



Figure 3. The parallel mean free paths from the Hussein 
& Shalchi simulations. Shown are the results obtained 
for slab/2D turbulence (dotted line), isotropic turbulence 
(solid line), Goldreich-Sridhar turbulence (dashed line), 
and NRMHD turbulence (dash-dotted line). 
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Figure 4. The perpendicular mean free paths from 
the Hussein & Shalchi simulations. Shown are the 
results obtained for slab/2D turbulence (dotted line), 
isotropic turbulence (solid line), Goldreich-Sridhar tur¬ 
bulence (dashed line), and NRMHD turbulence (dash- 
dotted line). 



Figure 5. The ratio Ax/Ay from the Hussein & 
Shalchi simulations. Shown are the results obtaiued for 
slab/2D turbuleuce (dotted line), isotropic turbulence 
(solid line), Goldreich-Sridhar turbulence (dashed line), 
and NRMHD turbulence (dash-dotted line). 



Figure 6. The perpendicular mean free paths \_l/L 
versus the parallel mean free paths A||/L from the Hussein 
& Shalchi simulations. Shown are the results obtained 
for slab/2D turbulence (dotted line), isotropic turbulence 
(solid line), Goldreich-Sridhar turbulence (dashed line), 
and NRMHD turbulence (dash-dotted line). 
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Figure 7. The perpendicular mean free paths Ax/^H 
versus the parallel mean free paths A||//|| for noisy slab 
turbulence. Shown are the simulations obtained by em¬ 
ploying the Hussein & Shalchi code (diamonds) and the 
analytical results (dotted line) represented by Eq. (29). 
For the sake of comparison we have also shown the sim¬ 
ulations for slab/2D turbulence (crosses). 
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Figure 8. The ratio Ax/Ay versus the parallel mean 
free paths A||/Z|| for noisy slab turbulence. Shown are the 
simulations obtained by employing the Hussein & Shalchi 
code (diamonds) and the analytical results (dotted line) 
represented by Eq. (28). For the sake of comparison we 
have also shown the simulations for slab/2D turbulence 
(crosses). 



Figure 9. The mean free paths for slab/2D turbulence. 
Shown are the parallel mean free paths from the Hussein 
& Shalchi simulations (dotted line) and PADIAN simu¬ 
lations (solid line) as well as the perpendicular mean free 
paths from the Hussein & Shalchi simulations (dashed 
line) and PADIAN simulations (dash-dotted line). 
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Figure 10. The mean free paths for isotropic turbulence. 
Shown are the parallel mean free paths from the Hussein 
& Shalchi simulations (dotted line) and PADIAN simu¬ 
lations (solid line) as well as the perpendicular mean free 
paths from the Hussein & Shalchi simulations (dashed 
line) and PADIAN simulations (dash-dotted line). 



Figure 11. The mean free paths for Goldreich-Sridhar 
turbulence. Shown are the parallel mean free paths from 
the Hussein & Shalchi simulations (dotted line) and PA¬ 
DIAN simulations (solid line) as well as the perpendic¬ 
ular mean free paths from the Hussein & Shalchi sim¬ 
ulations (dashed line) and PADIAN simulations (dash- 
dotted line). 



Figure 12. The mean free paths for NRMHD tur¬ 
bulence. Shown are the parallel mean free paths from 
the Hussein & Shalchi simulations (dotted line) and PA¬ 
DIAN simulations (solid line) as well as the perpendic¬ 
ular mean free paths from the Hussein & Shalchi sim¬ 
ulations (dashed line) and PADIAN simulations (dash- 
dotted line). 











